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ABSTRACT 

Abell 2218 is one of a handful of clusters in which X-ray and lensing analyses of the 
cluster mass are in strong disagreement. It is also a system for which X-ray data 
and radio measurements of the Sunyaev-Zel'dovich decrement have been combined 
in an attempt to constrain the Hubble constant. However, in the absence of reliable 
information on the temperature structure of the intracluster gas, most analyses have 
been carried out under the assumption of isothermality. We combine X-ray data from 
the ROSAT PSPC and the ASCA GIS instruments, enabling us to fit non-isothermal 
models, and investigate the impact that this has on the X-ray derived mass and the 
predicted Sunyaev-Zel'dovich effect. 

We find that a strongly non-isothermal model for the intracluster gas, which im- 
plies a central cusp in the cluster mass distribution, is consistent with the available 
X-ray data and compatible with the lensing results. At r< 1', there is strong evi- 
dence to suggest that the cluster departs from a simple relaxed model. We analyse the 
dynamics of the galaxies and find that the central galaxy velocity dispersion is too 
high to allow a physical solution for the galaxy orbits. The quality of the radio and 
X-ray data do not at present allow very restrictive constraints to be placed on Hq. 
It is apparent that earlier analyses have under-estimated the uncertainties involved. 
However, values greater than 50 km s~^ Mpc~^ are preferred when lensing constraints 
are taken into account. 

Key words: galaxies: clusters: individual: A2218 - X-rays: galaxies - dark matter - 
cosmic microwave background - gravitational lensing - distance scale 



1 INTRODUCTION 

The masses of galaxy clusters can be determined in three 
main ways: using the velocity dispersion of the galaxies, the 
pressure gradient in the hot intracluster gas derived from 
X-ray imaging and spectroscopy, and by analysing the lens- 
ing of background galaxies by the cluster potential. Each of 
these approaches involve assumptions and are vulnerable to 
various systematic errors. It is therefore useful to compare 
the results of the different techniques. The X-ray and lens- 
ing approaches are generally considered the most reliable. 
Results fro m them have now bee n compared for a number 
of clusters (Fort & Mellier 1994). The agreement is often 
reasonable, but there are a few spectacular exceptions, of 
which Abell 2218 (hereafter A2218) is the most weU studied 
example. 

A2218 is an optically compact (core radius « 1'; 
Dressi er 1978) cluster of galaxies, located at a redshift of 
0.171 ( [Kristian, Sandage fc Westphal 1978), and c lassified as 



metric studies (Pello-Descayre et al. 1988; Pello et al. 1995) 
suggest the existence of a second, smaller galaxy concen- 
t ration, displaced from the cD by 6 7". Spectroscopic study 
(Le Borgne, Pello & Sanahuja 1992), performed on the cen- 
tral region (< 4') of A2218, has provided redshift informa- 
tion for 66 of the objects within the core and shown that 
the average velocity dispersion is 1370 km s~^. 

A succession of X-ray telescopes have allowed the prop- 
erties of the hot gas within A2218 to be established. With 
the Einstein IP C fc HRI (jPerrenod fc Henry 198l[ ) and 
ROSAT PSPC dSiddiqui 1995| ), the emission was found to 
be smooth (on scales ~ 1 ), azimuthally symmetric and cen- 
tred on the cD galaxy. Fitting a polar profile of the surface 
brightness with a Kin g model gave a core radius of 58" lie 
and a /3- value of 0.63 (Boynton et al. 19821). Integrated spec- 



tral analyses with Ginga gave a gas temperature of 6.72 



+0.5 



keV and metallicity of 0.2I° j Z© (McHardy et al. 1990) 



richness class 4 (Abell, Corwin Jr & Olowin 1989). The clus- 
ter appears well relaxed, with the majority of the galaxies 
centred around the sole cD galaxy. However, detailed photo- 



By virtue of Ginga's bandwidth, this determination is com- 
monly accepted as the most accurate estimate of the mean 
gas temperat ure. Most recently , deep observations with the 
ROSAT HRI ( [Markevitch 1997 ) have shown the presence of 
significant X-ray substructure within the cluster core, sug- 
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gesting that the cluster may have undergone a recent merger 
event. This may account f or the absence of any sign s of a 
cooling flow in the cluster (Arnaud 1991; White 1996). 

Several previous com parisons ([Mirafda^Escude fc Babul 
1995|rR ' 
Kneib ] 



_b et al. 1995| [Kneib et al. 1996| ; [Natarajan & 



] 996| ) between strong lensing and X-ray analyses have 
found a factor of 2 discrepancy in the gravitating masses 
predicted by the two methods. Suggested explanations have 
centred upon the assumption of hydrostatic equilibrium for 
the cluster gas, the possibility that magnetic fields may pro- 
vide significant pressure support to the gas and the presence 

of subs truoturc within the oluatcr. 

Thb — Suujraev Zel'dovich — decromo nt associated with 



A2218 has been extensively studied (Jones et al. 1993 



Birkinshaw fc Hughes 1994 ; Saunders 1996 ). These results 
have been used, in coniunct ion with X-ray dat a, to con- 



strain l ihe lluhble cousCaul fcilk &. While 197a| , [ tvlcllaidj' 



et al. l t|90 | ; [Birkinohaw fe Hughoo 10oi[ )7 Thopo analyaoD hav 



been made in the absence of reliable information about tem- 
perature variations in the intracluster gas and have there- 
fore been forced to m ake simplifying assurnptions, such as 
that of isoth ermality (Mclfardy et al. 199(^ Birkinshaw & 
Hughes 1994; Kneib et al. 1995). This assumption is without 



a strong theoretical foundation an d conflicts with the results 
of most cosmological simulations ( NavarrOj_^Vonk & White 



1995 



^ avarro, Frenk & White 1996; Tormcn, Bouchet & 



White 



99(^ ) , which show temperature declining with radius, 
and mass distributions which have a central cusp. The ques- 
tion then arises as to whether simplifying assumptions have 
significantly biased the conclusions of previous X-ray anal- 
yses. For example, is the apparent discrepancy between the 
X-ray and lensing masses unavoidable or does it arise sim- 
ply from the use of inappropriate assumptions in the X-ray 
analysis? 

Motivated by the desire to avoid such restrictive as- 
sumptions, we have carried out an X-ray analysis which 
combines the capabilities of ROSAT and ASCA. The lim- 
ited spectral bandwidth and resolution of the ROSAT PSPC 
is compensated for by the superior spectral properties of 
ASCA. Conversely, the poor spatial performance of ASCA 
is complemented by the higher spatial resolution of ROSAT. 
This approach has never before been applied to A2218. 

The central aim of this paper is to improve our under- 
standing of A2218 by comparing the results of our X-ray 
analysis with lensing, SZ and galaxy velocity studies. It also 
serves as a case study on the possible dangers of assum- 
ing an isothermal gas, when one has no information to the 
contrary. Throughout the paper we assume an Einstein-de 
Sitter cosmology with Q=l, qo=0.5 and Ho — 50 km s~^ 
Mpc~^, except where otherwise stated. 



2 X-RAY ANALYSIS 

The objective of the analysis is to use spatially and spec- 
trally resolved X-ray data to constrain models of the distri- 
bution of gas properties in the cluster. For an in-de pth dis- 
cussion of the procedures covered in this section, see Cannon 
(1997)1 



2.1 Spectral-image modelling 

We work with X-ray spectral images, which constitute 
blurred records of the spectral properties of the cluster pro- 
jected along the line of sight. Since information about the 
disposition of material perpendicular to the plane of the sky 
is not available, it is necessary to make some assumption 
about the geometry of the source. We assume that the clus- 
ter is spherically symmetric. In prac tice, A2218 is slightly 
elliptical, with an axis ratio of 0.8 (Siddiqui 1995). How- 
ever this modest ellipticity sho uld not introduce anv serious 



errors into our derived masses ( Fabricant, Rybicki fc Goren- 
stein 198^ ). 



It is important to allow for the spatial and spectral 
blurring introduced by the telescope, as described by the 
instrument point spread function (psf) and ener gy response 



matr i x. We adopt a fo rward fitting approach ( [Eyles et al 



1991; Watt et al. 1992), in which the properties of the gas 



are parameterised as analytical functions of cluster radius. 
The emission from each spherical shell in the cluster is com- 
puted using a Raymond & Smith (1977), hereafter RS, hot 



plasma code. After correcting for cluster redshift, the spec- 
tral emissivity profiles are folded through the instrument 
spectral response, projected along the line-of-sight, rebinned 
into an xy grid and blurred with the psf. This produces a 
predicted spectral image which can be directly compared to 
the observed data, using a maximum-likelihood statistic. It- 
eratively adjusting the model parameters results in a best-fit 
to the data. 

Using analytical forms for the radial distribution of gas 
properties has the advantage of regularising the solution (i.e. 
suppressing instabilities in the deprojection and deblurring 
processes), however one runs the risk that the solution may 
be dictated by the mathematical function imposed. This can 
lead to overconfidence in derived results, as acceptable al- 
ternatives which might fit the data are ruled out by the lim- 
itations of the available models. The commonly employed 
restriction of isothermality is an extreme example of this. 
We attempt to avoid this problem by using a range of radial 
functions. This is particularly important for the tempera- 
ture and we use not only a number of parametric forms for 
Tgas{r), but also an alternative approach in which Tgas{r) is 
determined indirectly, by fitting a model for the mass distri- 
bution, as discussed below. The gas density profile is much 
more readily determined by the X-ray data, so we have fitted 
only two radial forms. 

Assuming that the intracluster gas is in hydrostatic 
equilibrium in the potential well of the cluster, the total 
gravitating mass within radius r, from the centre of the clus- 
ter, is related to the gas temperature and density by: 



Mgrav{r) = - 



kTgas ('^) 

Gfimp 



d\npgas{r) d\nTgas{r) 



din r 



dlnr 



(1) 



where Pgas{r) is the gas density, Tgasir) the gas tempera- 
ture, fj. the mean molecular weight and nip the proton mass. 



2.1.1 Gas density 

The gas density is well constrained by the X-ray surface 
brightness, since the PSPC is largely insensitive to variations 
in temperature for T> 3 keV. Surface brightness profiles are 
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generally w ell fitted by core- index type models (King 1962 
King 1972|): 



pgas(r) = pgasfi 1 + (r/r^) 



(2) 



where pgasfi is the central gas density normalisation (amu 
cm~^), Tc the core radius (arcmin) and Ofp the density in- 
dex (unitless). The main deviations from this form occur at 
small radii, where cooling flows give rise to surface bright- 
ness cusps in many clusters, though not in A2218. 

Recent N-body studies ( Navarro , Frcnk fc White 1995t 



Navarro, Frenk fc White 1996| ; [Tormen, Bouchet fc White 
199( ^) hj ave achieved good fits to dark matter (DM) and gas 
profiles in simulated clusters with an alternative descrip- 
tion. The profiles are found to steepen progressively, from 
Pgas{r) oc in the core, to near the virial radius, 
following the form 



Pgas{r) = pgasfi \x(l + x) 



(3) 



where x = r/vs, Vs being the scale radius (arcmin). Both 
of the above analytical forms have been fitted to the X-ray 
data for A2218. 



2.1.2 Gas temperature 

The gas temperature distribution is less well determined, 
since this requires a combination of spatial and spectral res- 
olution which has not generally been available in the past. 
We consider a variety of simple models: a linear temperature 
ramp (LTF), 



Tgas{r) = Tgasfi — Pr 



(4) 



where Tgasfi is the gas temperature (keV) at the cluster 
centre, /3 the temperature gradient (keV arcmin"^) and r 
the radius (arcmin); a King- type temperature description 
(KTF), 



21-/3 



(5) 



Tgas{r)= Tgasfi [l + [r/r^f] 

where Tt is the temperature core radius (arcmin) and /? the 
temperature index (unitless); and a polytropic temperature 
description (TTF), 



Tgas{r) = Tgasfi [l + (r/rc) 



21 ap(l-7) 



(6) 



where r^ is the gas density core radius (arcmin) and 7, the 
polytropic index (unitless), is fitted as a free parameter vary- 
ing between isothermality (7=1) and adiabaticity (7=5/3). 

2.1.3 Gravitating mass 

An alternative to fitting Pgas{r) and Tgasir) is to fit Pgas{r) 
and yigravir). The corresponding temperature profile can 
then be inferred, via Equation ^. We use several alternative 
f orms, motivated by the distribution of galaxies in clusters 
(Rood et al. 1972), and by the results of N-body simulations. 
These include: a core-index description (DMF), 



pDM(r) = Pdm.o [1 + [r/rcf 



(7) 



where Pdm.o is the central dark matter density normalisation 



(amu cm 



the core radius (arcmin) and qdm the den- 



sity index (unitless) ; a model ba sed upon the simulations of 
Navarro, Frenk fc White (1995)| (DNF), 



PDM(r) = Pdm.o [x{l + a;)^] ^ (8) 

where x = r/rs and is the scale radius (arcmin); and a 
Hernquist profile (DHF), 

PDM(r) = Pdm,o/ [27rfe(l + bf] (9) 
where b — r/rs and rs is the scale radius (arcmin). 

2.1.4 Fitting the models 

Determination of the best-fit parameters for a cluster model 
proceeds in the way commonly employed for spectral fitting. 
The fit statistic employed is maximum likelihood, rather 
than chi-squared, since the data are generally strongly Pois- 
sonian. The fit and its local slope are determined at some 
initial position in the parameter space. This information is 
used to predict an improved set of model parameters and 
the fit statistic re-determined. The process is iteratively re- 
peated until the statistic slope falls below a pre-determined 
value. One limitation of this method is, however, that the 
fitting tends to follow the local gradient in the statistic, until 
it encounters a minimum. Thus the fit can become trapped 
in a "valley", which it regards as the best-fit result, even 
though a more suitable combination of parameters may oc- 
cur elsewhere. To avoid this, we randomly perturb models 
during analysis and force them to re-fit (to check if the same 
minimum is produced). 

Confidence regions can be derived for each best-fit pa- 
rameter, by offsetting the parameter of interest from its best- 
fit value (both above and below the best-fit), and reoptimis- 
ing the other parameters. The resulting increase in the fit 
statistic, from its optimum value, is used to determine what 
offset would need to be applied in order to create a user- 
defined change in the statistic. This defines the required 
confidence interval. We use the form of the maximum likeli- 



hood statistic introduced by Cash (1979), such that changes 



in the statistic have the same significance as changes in chi- 
squared. Hence, for each parameter, an increase in the Cash 
statistic of 1 corresponds to a 68% confidence interval, and 
an increase of 2.71 to 90% confidence. The above process is 
repeated for each model parameter for which an error esti- 
mate is required. 

Errors in physical quantities which are functions of ra- 
dius (such as mass or temperature) are generally affected 
by several model parameters. We derive error envelopes for 
such quantities by taking the outer envelope of all of the 
curves generated by perturbing each free parameter to its 
upper and lower error bounds. Because these envelopes are 
derived using every parameter combination, each offset to 
their error bounds, the result is a conservative estimate of 
the statistical uncertainty. 

Once the total gravitating mass distribution has been 
determined (using Equation |l| if the mass has not been mod- 
elled directly) the various mass components in the cluster 
can be separated. The gas mass profile is calculated from 
the fitted parameters. The galaxy mass profile can be con- 
structed from the observed luminosity profile for the cluster, 
assuming a constant mass-to-light ratio. Subtraction of these 
components from the total mass profile then yields the dark 
matter profile. 
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Figure 1. Gas temperatures derived in projected annuli. The 
barred crosses are taken from ROSAT PSPC analysis and the di- 
amonds from ASCA GIS analysis. The results suggest a temper- 
ature decline with radius, with weak evidence (from the PSPC) 
for central cooling. 



a full spatial and spectral analysis, the data were formed 
into images of channel width 10, over the channel range 11- 
230 (approximately 0.2-2.3keV). This results in a data cube, 
from which specific regions can be selected and analysed. 

Within the cluster emission there is one bright source 
contaminating the data, located at a radius of 11.1' from 
the cluster centre. In the PSPC analysis, the point source 
can be eliminated by ignoring the data collected in that 
region. However, this is not possible with ASCA, since the 
extended psf ensures that it is not discernable as a discrete 
source. To determine whether the source significantly affects 
our analysis, we model it using the PSPC data. The cluster 
emission is first fitted with the point source pixels removed. 
This model of the cluster emission is then subtracted from 
the original data, leaving behind just the point source, which 
is fitted with a power-law spectral model. 

The fitted index, 1.25, indicates the softness of the 
source - the majority of the fiux is emitted below 0.5keV. 
This is consistent with identification of the source as 



SA0 17151, a bright star with a soft spectrum ( Markevitch 



1997). If the PSPC-determined model for the point source is 
subtracted from the ASCA data, the cluster fits are modified 
to the extent that a 1.5% difference in the total gravitating 
mass at 2Mpc results. This effect is negligible compared to 
other errors, so no attempt has been made to remove the 
source from the ASCA data. 



2.2 ROSAT PSPC reduction 



2.3 ASCA GIS reduction 



The aim of the ROSAT analysis is to obtain well constrained 
gas density parameters, which can then be utilised in the 
ASCA analysis. 

The raw data, obtained on May 25th 1991, were re- 
duced using the Starlink ASTERIX X-ray analysis package. 
Periods of high background were removed from the data, 
reducing the effective exposure time to 42 ksec but making 
the background subtraction substantially more reliable. 

Subtraction of the X-ray background was accomplished 
by selecting data from an annulus (27' — 33'), ignoring pix- 
els covered by the detector support structure or contain- 
ing point source emission. This background sample was 
then extrapolated to cover the whole field, using the PSPC 
energy-dependent vignetting function. Since the X-ray sur- 
face brightness profile for A2218 can be traced to a maxi- 
mum radius of 12', the chosen background annulus is free of 
source emission. In the following analysis, the data are re- 
stricted to lie within 9' to avoid possible systematic efi'ects 
from uncertainties in the background subtraction at large 
radius. 

The exposure-corrected, background-subtracted PSPC 
data were summed to provide an integrated cluster spec- 
trum and split into concentric annuli, centred on the clus- 
ter core, within which spectra were extracted. Fitting these 
spectra gives an indication of the depth of the cluster poten- 
tial well and the radial structure of the cluster gas density 
and temperature parameters (see Fig 0). However, any gra- 
dients present in these distributions will tend to be underes- 
timated, due to the smoothing effects of the instrument psf 
and projection along the line-of-sight. Using annuli of width 
greater than the instrument psf minimizes the former effect. 

In order to create a spectral image dataset, allowing 



The ASCA analysis aims to constrain the gas temperature 
and metallicity profiles, using the PSPC derived gas density 
profile. A2218 was observed by the ASCA X-ray telescope 
on April 30th 1993. In this paper we use only data from 
the two gas imaging spectrometers (GIS2 and GIS3), since 
these have a wider field of view (~ 50' diameter) and greater 
high-energy detector efficiency (up to lOkeV) than the CCD 
detectors (SISO and SISl). An additional reason is that an 
accurate model for the large, asymmetrical and energy de- 
pendent psf is available for the GIS detectors, constructed 
from Cyg X-1 observations. These restrict analysis to a max- 



imum radius of 1 8' and an energy-range of 1.5-11 keV ( Taka-^ 
hashi ct al. 199^ ) , which is not a significant limitation in the 
case of A2218. 

Standard proced ures for ASCA ana lysis, followed in this 
paper, are given by Day et al. (1995). The recommended 



screening criteria are applied to the raw data, removing data 
taken during times of high background fiux. Subtraction of 
the X-ray background is complicated by the telescope psf. 
This has the effect of ensuring that no region of the detector 
is free from source fiux. Hence, we extract an "average" back- 
ground da taset from the pu blicly distributed set of blank-sky 
pointings (Day et al. 1995). These datasets suffer from mild 



point-source contamination, as sources are not completely 
averaged out, but are currently the best available solution. 

The results of a naive annular spectral analysis of the 
ASCA data are shown in Fig ^ However, it is important 
to bear in mind the limitations o f this approach. Cross - 
talk between annuli is significant ( |Takahashi et al. 1994 ), 
and the energy-dependent spreading of flux results in a dis- 
torted temperature profile, such that analysis of a simulated 
isothermal cluster would give a temperature which rises with 
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radius. In practice, the temperature appears to decline with 
radius, indicating that a real gradient is present. 

For 3D analysis, spectral-image datasets with contigu- 
ous energy bands of width 50 raw channels are created. Since 
the cluster centre is offset from the detector centre, data be- 
yond a radius of 9' are not fitted (the offset added to the 
radial extent of the source is similar to the maximum radius 
where psf calibrations apply) . This restriction minimises the 
effects of poor calibration and high background near the de- 
tector edge. As both GIS instruments behave similarly, the 
datasets are fitted simultaneously. The Cash statistic has 
been used to identify best-fit models, but similar results for 
best-fit parameters and for comparison between the quality 
of fit of different models, is obtained using the statistic. 



3 R ESULTS 

We first compare the results of our analysis with published 
studies, to investigate whether the fitted models are consis- 
tent with earlier work on A2218. 

Integrated spectral analyses of clusters produce "mean" 
quantities which are r epresentative of the enti re object. As- 
suming isothermality, McHardy et al. (1990) derived a gas 
temperature of 6.72^q'4 keV together with an iron abun- 
dance of 0.21q'2 Zq (using an RS emission code, where 
all other heavy element abundances were fixed at 0.5 Zq). 
This is in agreement with an earlier, much les s well con- 
strained examination (Perrenod fc Henry 1981). More re- 



lier, disrupted cooling fiow (as noted by |Siddiqui 199f^ ). The 
ASCA analysis suggests that the metallicity may be slightly 
lower than the integrated value in the cluster centre with a 
sha llow radial rise. Howev er, all points are consistent with 
the [McHardy et al. (1990)| value of 0.2 Zq. 

Analysis using spectral-image datasets allows extrac- 
tion of the 3D gas density and temperature distributions 
within A2218. In the subsequent cluster analysis, both 'tem- 
perature models' (fitting for Pgas{r) and Tgas{r)) and 'mass 
models' (fitting for Pgas{r) and Mtot{r)) are used. Parame- 
ters representing the metallicity and cluster position are also 
fitted. 

The best constrained parameters derived from PSPC 
analysis pertain to the shape of the gas density distribution. 
Comparable analyses have been carried out us i ng Einstein 
([Perrenod fc Henry 1981 ; Boynton et al. 1982; Birkinshaw 
Hughes 1994| ), and ROSAT ( [Siddiqui 1995| ) data. These 



studies agree that when a King profile is assumed, A2218 is 
well modeled with a core radius, rc, of ~ 1' and an index. 
Up, of 1 (equivalent to a /3-value of 0.67). Higher resolu- 
tion analysis, using th e R OSAT HRI, has bee n performed by 



Squires et al. (1996) 



cently, Mushotzky fc Loewenstein (1997)| have derived T = 
7.2 keV and Z = 0.18 ± 0.07 Zq, from an integrated ASCA 
spectrum. 

Fitting an RS model to our ROSAT data results in a 
gas temperature of 4.71q g keV and a hydrogen absorption 
column of 2.61q^x10^° cm~'^ (with metallicity fixed at the 
Ginga value). This abs orption agrees wit h the Stark level of 



2.58Io igxlO^" cm"-^ ( Ptark et al. 19921). The temperature 



is lower than the Ginga result of McHardy et al. (1990) 



but the energy range of the PSPC is not very suitable for 
determining the t emperature of such hot gas. It has been 
found previously ([Markevitch fc Vikhlinin 1997|) that PSPC 



and [Markevitch (1997)| . These studies 
detect the presence of smaller scale structure, with three sur- 
face brightness peaks visible within the central arcminute, 
none of which coincides with the central cD galaxy. This 
structure cannot be resolved with either the PSPC or GIS 
detectors, and indicates that the core of A2218 departs from 
a fully relaxed state. 

Since ASCA is poor at constraining the gas density dis- 
tribution for such a compact source, the PSPC fit values 
for core radius and index are carried over into the ASCA 
spectral-image analysis. To obtain the appropriate density 
parameters, a linear temperature ramp model is fitted to 
the entire PSPC data within a radius of 9 arcmin, with 
the temperature parameters fixed at those derived from an 
initial fit to the GIS data. This model is then re-fitted to 
the GIS data, complete to a radius of 9 arcmin, with the 
gas density core radius and index fixed, allowing a fit of 
the temperature parameters. In an iterative process, this 
model is alternately fitted to the PSPC and GIS datasets 
until no further change is observed in the parameter values. 
With consistency achieved, the "standard" gas distribution 

g;Q3 arcmin, ap — 



0.96: 



[, in good agreement with the comparable results 



discussed above. If an NFW profile is assumed, the required 



Mushotzky & Loewenstein (1997) 



results. 



results tend to be biased low for high temperature clusters. 
Fitting an integrated spectrum simultaneously to the GIS2 
and GIS3 data gives a gas temperature of 6.73j;Q 44keV and 
a metallicity of 0.20j]Q ggZ© (with the hydrogen column fixed 
at the PSPC value), in good agree ment with the Ginga and scale radius, r,, is 10.27+^:20 arcmin. However, the NFW 

parameterisation is neither preferred nor disallowed by the 
PSPC data, so we only use a King parameterisation for the 
gas density distribution in the following analysis. 

3.1 All cluster data 

Fixing the gas density shape parameters at these King val- 
ues, a range of models were fitted to the ASCA spectral- 
image data. On the basis of their Cash statistic, a set of 
models best-fitting the observed cluster data within a ra- 
dius of 9 arcmin is selected (Table ^ lists the main parame- 
ters for each model) as representative of A2218. As can be 
seen from the Table, the isothermal model is a significantly 
poorer fit to the data. It is included for comparison with 
the other models. All of the temperature profiles obtained 
for this set (apart from that for the isothermal model) are 



The ability to extract and analyse spectra from inde- 
pendent regions of the cluster represents a significant ad- 
vance over analysis of the integrated emission. Fig ^ shows 
the derived annular temperature profiles from both ROSAT 
and simultaneous GIS2/GIS3 anal ysis. The PSPC results 
are consistent with Siddiqui (1995), with a possible temper- 
ature drop visible in the central bin. However, the evidence 
for central cooling is statistically rather weak, and the de- 
rived central cooling time of ~ l.SxlO^^yr is comparable 
with the Hubble time, so a strong steady-state cooling flow 
appears to be ruled out. 

If the hydrogen column is fitted, using the PSPC data, 
it is found to be consistent with the Stark value dtark et al 



1992) throughout the cluster, apart from a slight rise in the 



centre. This may be due to matter deposited by an ear- 
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plotted in Fig |^ Beyond the central arcminute (^^230 kpc) 
the profiles are in good agreement and quite non-isothermal. 
The typical 90% confidence envelope for an individual model 
(remember that these are conservative envelopes) generally 
encompasses the spread of these best-fit profiles (a single 
error envelope is included to illustrate this point). Within 
the central region, which lies within a single ASCA psf, a 
greater spread in temperature is allowed by the data. How- 
ever, despite this divergence in temperature at small radius, 
all of the best-fit models, bar the isothermal model, have 
similar Cash statistics (see Table 0). 

The abundance of heavy elements has been assumed, 
in all of the above models, to be constant over the cluster. 
However, the spectral capabilities of ASCA allow us to test 
this assumption. Allowing a linear metallicity gradient with 
the best-fit linear temperature ramp model gives a slope of 
S.Slg'g X 10~^ Zq arcmin"^, a value consistent with uni- 
form metallicity. The effect of this best-fit slope upon other 
model parameters is negligible, hence freezing the metallic- 
ity gradient at zero does not bias our analysis. 



3.2 Central cluster data 

We now examine the claim that lensing analyses require a 
larger cluster mass than is consistent with the X-ray data. 
The model which provides the greatest gravitating mass at 
the critical radius (~85 kpc) is therefore selected for exam- 
ination. This model, the DHF best-fit, using a Hernquist 
description for the DM distribution and a King description 
for the gas density profile (see Table implies a high cen- 
tral gas temperature. In the following comparisons, we refer 
to this model as the "maximum- mass model" (MMM). 

The MMM temperature profile features a factor of 2 rise 
within the central 1'. This raises two questions: does such a 
steep temperature gradient raise physical problems (would 
it be convectively unstable?), and is it consistent with the 
X-ray spectral data observed within the central region? We 
will return to the first question in Section ^. 

To address the latter question, GIS3 spectra integrated 
within r = 1', where the MMM temperature rises steeply, 
are compared with the predictions of the MMM and isother- 
mal models in Fig ^ For clarity, only the energy range of 
5-10 keV is displayed, since this is where the impact of very 
hot gas will be most apparent. Although the MMM provides 
a reasonable fit to the spectral imaging data as a whole, it 
does not follow that the data in the central regions need 
be consistent with the high model temperature. In practice, 
for both instruments (while the GIS2 spectra are not shown 
in Fig c , they behave similarly to the G1S3 spectra) the 
MMM IS a good match to the data, t he isothermal model 
is, however, also consistent with the restricted dataset. The 
conclusion from this is that while the data do not rule out a 
central temperature rise, they do not require one either. The 
reason for this is that within the ASCA bandpass, plasmas 
with temperatures of 8 and ISkeV do not have substantially 
different spectral signatures and are thus difficult to differ- 
entiate (this is analogous to the difficulty that the PSPC has 
in dealing with clusters hotter than 2-3keV). This problem 
is compounded by the limited spatial resolution of ASCA. 



Model 


Central 




Relative 


Form 


gas 


Gradient 


Cash 




(keV) 


(various) 


Statistic 


LTF 


9.64 


0.63 


-6'^6 08 


LTF (ISO) 


7.90 


zero 


-626.85 


TTF 


19.01 


1.20 


-634.53 


KTF 


8.98 


3.18 


-635.96 


DNF 


15.59 


n/a 


-636.12 


DMF 


9.24 


n/a 


-635.20 


DHF (MMM) 


20.42 


n/a 


-635.33 


MMMC 


16.44 


n/a 


-632.50 



Table 1. Temperature and statistic parameters for a variety of 
fitted 3D models are shown. Note that the MMMC is not a best- 
fit model in the same sense as the others, as is described in the 
text. For all models, the gas density core radius, Tc, and index. 
Op, are fixed at the values determined from the ROSAT data. 
Explanations for the model acronyms are given in the text. 



4 COMPARISON OF X-RAY AND LENSING 
MASSES 

Deep optical observations of A2218 reveal a number of major 
arcs and a wealth of minor arclets. These features are the 
result of gravitational lensing, an effect which is independent 
of the physical state of the cluster gas. Instead, uncertainty 
lies in the characteristics of the background galaxies and the 
possibility of matter sub-clumps along the line-of-sight. 

Using the models fitted to A2218, we can derive pro- 
jected gravitating mass profiles, suitable for comparison with 
the results of lensing analysis. This involves assuming a max- 
imum outer radius for the cluster mass distribution. In the 
following analysis we take this to be the maximum radius of 
the data used in cluster fitting, which is 9' (~2 Mpc). The 
choice of projection radius has a minor impact upon the 
derived projected mass profiles (compared to other uncer- 
tainties) so long as the chosen radius is sufficiently large, 
> 2Mpc. For example, the difference in projected mass 
within 2Mpc, between models with maximum radii of 2.0 
and 3.0 Mpc, is 9%. 

Since lensing analysis measures the matter distribution 
on both small (strong lensing occurs where the surface mass 
density is high) and large (weak lensing is theoretically ob- 
servable to the edge of the cluster) scales, comparison with 
X-ray results is extremely informative. In Fig ^ the X-ray 
derived projected mass profiles for the MMM and isother- 
mal model are plot t ed, together with the len sing results of 



(1996) 



Kneib et al. (1995), Loeb fc Mao (1994) and Squires et al 



4.1 Strong lensing 

Strong lensing analysis has been performed, using spectro- 
scopically observed arcs, by a variety of groups (|Loeb & 



Mao 1994|; |Kneib et al. 1995|; [Kneib et al. 1996|; garaniti 



Petrosian & Lynds 199(; [Na^arajaiTTn^neib 1996 ). These 



analyses provide the gravitating mass within the critical ra- 
dius associated with the formation of giant arcs, which in 
A2218 is 22.1" (~ 85kpc). 

The two strong lensing points, plotted in Fig ^ at this 
radius, are consistent with a projected mass of ~6xlO^^M0, 
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Figure 2. Fitted 3D temperature profiles for the cluster gas. The MMM (solid) is plotted together with its 90% error envelope (dash- 
3dot) and the remaining best-fit models (dotted). All of the profiles feature gas temperatures which decrease by a factor of ~ 2 or more 
within the region of analysis, and are consistent beyond the cluster core. 



des pite their use of different mass distributions (|Kneib et al 



1995 _as jume a bipolar mass model while Loeb fc Mao 1994 et al. (1995) 



use an isothermal sphere distribution). The corresponding 
projected mass for the MMM is ^.'itW^lQ^^ Mq. The un- 
certainty associated with this mass is large enough to in- 
clude the lensing derived value, so, within the calculated 
errors, the X-ray and strong lensing analyses are consistent, 
when the MMM is used. The corresponding mass from the 
isothermal model is 2.81q'2^10^^ Mq (see Fig ^), in agree- 
ment wit h the factor of 2 discrepancy reported by previous 
analyses ([Loeb fc Mao 1994 [Miralda-Escude fc Babul 1995 
Kneib et al. 1995D 



The bipolar mass model of Kneib et al. (1995) 



was used 

by these authors to extract a further two masses (see Fig |^ , 
at the radius of the second mass clump (256 kpc) and at 
the maximum radius where giant arc constraints can be ap- 
plied (383 kpc). The projected, enclosed gravitating mass 
at these radii for both the MMM and the isothermal model 
are given in Table |^. No error estimates are given by Kneit 
et al. (; 995) for the lensing masses. However, it can be seen 



that both of these outer lensing points lie just outside the 
statistical error envelope of the MMM, unlike the mass de- 
rived at the critical radius. This is due to the MMM profile 
flattening at large radius. Note however that both of the 
outer predicted ma sses depend upon the bimodal model of 
Kneib et al. (1995), whereas the inner point can be derived 



in a model-independent fashion. 



The Kneib et al. (1995) bimodal mass model has re- 
cei ved additional suppo rt from detailed optical observations 
by Ebbels et al. (1996). In this analysis, the redshift of one 
of the faint arcs was determined for the first time and found 



to be in good agreement with the value predicted by Kneib 



The conclusion which can be drawn from this X- 
ray /strong lensing comparison is that it is not clear that 
the results from the two approaches are inconsistent. We 
have shown that when a model such as the MMM is used, 
which includes a central mass cusp, the predicted masses at 
the critical radius agree within the X-ray statistical error 
envelope (see Table H). If the gas is assumed to be isother- 
mal, the previously noted discrepancies can be reproduced 
(see Fig ^). 



4.2 Weak lensing 

A statistical an alysis of the weakly l ensed arclets has been 
carried out by Squires et al. (1996), allowing the slope of 
the surface mass density to be mapped within a radius of 
~ IMpc. Fig ^ shows that, within errors, the magnitudes 
of the weak lensing points are consistent with the MMM. 
However, the trend indicated by these points is for a steeper 
profile, providing more gravitating mass at large radii than 
predicted by the MMM. 

A sig nificant difficulty with the weak lensing analysis. 



noted by Squires et al. (1996), is the procedure used for 



normalisation. This is done by defining a reference annulus 
at large radius within which the cluster mass contrib ution is 
assumed to be negligible. Since the annulus used by Squires 



et al. (1996) has an inner radius of 800 kpc, while the X-ray 



surface brightness profile can be traced to r > 2 Mpc, some 
shear signal from the cluster must actually be present within 
the reference annulus. Hence the recovered normalisation is 
an underestimate, such that the weak lensing mass profile 
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Source/Mass 85kpc (xlO^M©) 256kpc (xlO^^M©) 383kpc (xlO^^M©) 

Kneib 1995 6.1 2.7 4.5 

Loeb&Mao 1994 6.4 

MMM 5.3l^;5 2.0l?;^ 3.0l°-J 

Isothermal S.stgj l-6lo? ^.S+qI 



Table 2. Projected gravitating masses from strong lensing and X-ray analysis are compared at the radii used by Kneib et al. (1995) 



At 



the critical radius the isothermal mass is a factor of 2 too low, although the discrepancy reduces greatly at larger radii. 




Energy (keV) 

Figure 3. The GIS3 spectrum from the central 1' (crosses) is 
compared to the prediction from the MMM model (solid) and the 
isothermal model (dashed). The MMM model provides more flux 
at high energy, because of the centrally increasing temperature, 
but the difference is very small. Both models are consistent with 
the observed data. 



can only be regarded as a lower bound to the projected mass 
The impact of this is examined further in S ection 



correction for this effect has been estimated by Squires et al 
(1996)[}o be a factor of ~1.2-1.6 in projected mass 

If the weak lensing points are adjusted to take ac- 
count of this factor, this has two important implications for 
the MMM comparison. First, the normalisation of all data 
points increases, bringing the inner points into better agree- 
ment with the MMM profile while the outer points move 
further from consistency. Second, because this adjustment 
is a DC effect for the reconstructed projected surface mass 
density (of the cluster), it does not act equally on the ra- 
dially integrated points. Hence, when the normalisation is 
raised, the slope of the weak lensing mass profile in Fig ^ 
increases, bringing it into greater conflict with the MMM 
profile. Thus the flatter slope of the MMM cannot be made 
consistent with the weak lensing points by shifting the ref- 
erence annulus to larger radii. 

Even when the normalisation adjustment is applied, this 
is insufficient to e nsure full consistency with the outer strong 



dieted at the same radii by weak lensing. There is, however, 
a well understood effect whereby the weak lensing signal is 
suppressed at small radii due to contamination by cluster 
galaxies, red u cing the derived w eak lensing mass (Kaiser & 



Squires 1995; Squires et al. 1996). Correcting for this would 



bring the inner weak lensing points points into greater con- 
sistency with both the strong lensing and X-ray results. 

Overall, then, this comparison indicates that at large 
radii, > 250 kpc, the weak lensing and X-ray analyses are 
reasonably consistent, though the weak lensing results tend 
to give more mass at large radii. Further work (observations 
of increased numbers of arclets over a wider field) is required 
to reduce the uncertainty in the weak lensing analysis and to 
move the reference annulus to larger radii, where the cluster 
mass contribution is lower. 



5 GAS ENTROPY 

In Fig P the derived gas entropy profiles for the best-fit clus- 
ter models are plotted. These indicate that beyond a radius 
of ~ 120kpc the entropy increases with radius, as is ex- 
pected for gas which is convectively stable. However, within 
this radius three of the best-fit models (including the MMM) 
exhibit a slight inward rise in entropy, making the gas con- 
vectively unstable. For the MMM, the predicted increase in 
entropy which occurs between a radius of 120kpc and the 
cluster centre is ~ 30%. The reason for this behaviour is 
that the temperature increases rapidly at small radii while 
the gas density is forced to flatten (due to the use of a King 
model description). However, it is precisely this temperature 
increase that allows the model to achieve consistency with 
the gravitational mass derived from strong lensing. 

As noted in Section ^ the NFW gas parameterisation 
is neither preferred nor disallowed by the PSPC data. How- 
ever, if the King gas density parameterisation used for the 
MMM is replaced by the NFW parameterisation, the central 
entropy drops by ~ 75% (See Fig ^). This difference is large 
enough to ensure that entropy rises continuously at all radii, 
removing the problem of convective instability. 



6 CENTRAL GALAXY MASS 

An alternative method for bringing the strong lensin g and X- 
ray ana lyses into agreement has been examined by 



Makino 



lensing points of Kneib et al. (1995) 

lie above the values of 1.5x10^'' Mq and 2.5x10'* Mq pre- 



which (see Table 



(1996), In this study, a massive central galaxy was embed- 
ded within the cluster to increase the predicted core mass 
without violating observational X-ray constraints. As the 
cD galaxy envelope can be traced to at least 25" (96kpc), 
which encompasses the critical radius, this is potentially an 
important consideration. 
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Figure 4. Projected gravitating mass profiles for tlie MMM (solid line) and isother mal (dotted line) ar e plo tted. Also shown is t he 90% 



error envelope for the MMM (dashed lin csl. Overlaid are stro ng lensing points from Kneib et al. (1995) and Loeb & Mao (1994) (boxes) 
together with the weak lensing points of ^quires et al. (1996) (crosses). 




Figure 5. Derived entropy profiles for the cluster gas. The MMM is plotted for both a King (solid) and an NFW (dashed) gas density 
distribution. The remaining best-fit models (dotted) are plotted using the King parameterisation. It can be seen that all of the models 
are consistent with a gas which is convectively stable, except within ~ 120kpc, where the high central gas temperatures of certain models 
(such as the MMM, see Table |l|) lead to an increase in entropy. Where the NFW gas density form is used, the derived entropy increases 
at all radii. 
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Makino (1996) tested this hypothesis by, firstly, con- 
structing a total mass distribution which included both clus- 
ter and cD components. These were parameterised by a King 
mass model (equation ^ with qdm = 3/2) and an isothermal 
sphere mass model (odm = 1), respectively. Secondly, the 
gas temperature distribution was constrained to be isother- 
mal at large radius (beyond the cD galaxy) and to rise or fall 
linearly with r at small radius. By assuming the gas to be in 
hydrostatic equilibrium, the gas density distribution corre- 
sponding to diflFerent temperature models was extracted and 
compared to the observed Einstein surface-brightness data. 
Constr aining the temperature t o be isotherma l at r > 



31", at the McHardy et al. (1990) 



derived value. 



Makinc 



(1996)[]found that models consistent with both the X-ray 
surface brightness and with enough central mass to account 
for the strongly lensed arcs had two notable features. Firstly, 
a cD component was required to ensure that the critical sur- 
face mass density was attained by the model. Secondly, the 
gas temperature rose sharply within 31", reaching a central 
value of ~ llkeV. Models which did not include a cD mass 
contribution, or which assumed the gas to be isothermal at 
all radii, were unable to provide sufficient mass at small radii 
to account for the lensing data, whilst remaining consistent 
with the Einstein data. 

The projected gravitating mass within 1.5 Mpc, from 



the model favoured by [Makino (1996)|, of 1.3x10^^ M© is in 



reasonable agreement with the MMM value of 9.0_5'3xl0 
Mq. This c onsistency is al so found at smaller radii, r=100 



kpc, where Makino (1996) 



compared to Q.G^^^^IO^ 



obtains a mass of 5.2x10 M© 
M0 for the MMM. It should be 
noted that the MMM employs a mass profile with a central 
mass cusp (pdm (x r~^) but does not include any discrete 
component corresponding to the central galaxy. We now ex- 
plore the effect of adding such an additional central compo- 
nent to our model. 

Adding a central galaxy of radius lOOkpc and adjustable 
mass normalisation to the MMM, we find that the statistical 
quality of the fit deteriorates. The largest additional mass 
allowed, at the 95% level, is 1.7 x 10^^ M© (within lOOkpc). 
From here on this model is referred to as the MMMC, since 
it represents the MMM plus a central galaxy mass. Note 
however that the MMM, which contains no central mass 
component, is still statistically preferred (see Table p]). 



This central mass is substantially less than the Makinc 
(1996)[]/alue of 5.2x10"Mq (recall, however, that our clus- 
ter profile contains a central cusp, whilst Makino's has a flat 
core). Compared to the MMM, it provides a mass increase 
at the critical radius of only ~5%. However, whilst the addi- 
tion of a central mass component has little impact at small 
radius, it has the rather counter-intuitive effect of providing 
more mass at large radius. The reason for this is that, in the 
case of the MMM, the DM distribution has a scale radius of 
1.89I"- jg arcmin. When a central mass component is added, 
the DM profile no longer needs to peak so sharply and its 



+0.95 
-1.49 



scale radius increases to 6.79 

Fig ^ shows the MMM and MMMC mass profiles com- 
pared to the strong and weak lensing results. It can be seen 
that the increased mass of the MMMC at large radius (com- 
pared to the MMM) allows it to be consistent, within its 
90% statistical co nfidence envelope, w ith both the two outer 
points de rived by iKneib et al. (1995) and the weak lensing 



points of Squires et al. (1996) 



We can use our models to estimate the baseline error 
involved in the weak lensing analysis, as a result of clus- 
ter mass residing wi thin the reference annulus employed by 
Squires et al. (1996)| . In Fig | the result of using the MMMC 



to correct the weak lensing results for this mass is also plot- 
ted. It can be seen that the points come into better agree- 
ment with the MMMC profile, particularly at r>500 kpc. At 
smaller radii the weak lensing points still lie systematically 
below the MMMC profile, although this is to be expected 
(see Section [l.2| ) due to dilution of the lensing signal by 
cluster galaxies. 

In Fig 1^ the derived gas entropy profile for the MMMC 
is compared with those from the best-fit X-ray models within 
the central 500kpc region, where the MMM shows a notice- 
able rise. It can be seen that the MMMC provides a flatter 
central entropy distribution, such that the gas is less likely 
to be subject to convective instability. 



7 ANALYSIS OF GALAXY MOTIONS 

The dynamics of cluster galaxies provide a further way of 
investigating the mass distribution in clusters. In the case of 
A2218, galaxy redshifts are available only within a few core 
radii of the centre. How ever this includes th e region where 

and the high 



the lensing analysis of Kneib et al. (1995 



resolution X-ray observations of Markevitch (1997), suggest 



that the potential may be seriously disturbed. 

Galaxy velocity and position data have been obtained 
from the NASA Extragalact ic Database. This consists 
largely of the data obtained by Le Borgne, Pello fc Sanahuja) 



(1992), who performed an extensive photometric survey of 
th e cluster core. On the basis of the 3(j clipping technique 
of [Vahil fc Vidal (1977)1 49 of the 53 galaxies with mea- 
sured redshifts are identified as cluster members. All of these 
galaxies lie within 3' of the X-ray centroid of the cluster, 
which has been adopted as the centre for this optical analy- 
sis. The average line of sight velocity disper sion, calculated 



using the 49 cluster galaxies, and corrected ( Harrison 1974 ) 



to the cluster rest frame, is 1354_jj^g km s~ . 

This galaxy distrib ution has been studied usi ng the 

using 



techniques described by Hobbs & Willmore (1997) 



the Jeans equation to relate the spatial and velocity dis- 
tribution of the galaxies to the gravitating mass profile for 
the cluster. The aim of such an analysis is to determine 
the radial behaviour of the anisotropy in the galaxy veloc- 
ity distribution, since no information about galaxy orbits is 
otherwise available. Since the radial distributions of galaxy 
density and velocity are projected along the line of sight, the 
analysis requires an extrapolation to large radius of the X- 
ray determined mass proflle, the galaxy velocity dispersion 
profile and the galaxy surface density profile. Extrapolating 
the galaxy profiles is particularly uncertain because data 
only extend out to ~ 2' (~ 500kpc) and include no informa- 
tion (except in projection) regarding the behaviour of these 
profiles beyond this region. 

The anisotropy is studied through /3, the anisotropy pa- 
rameter, which is 1 in the case of purely radial orbits, for 
an isotropic velocity distribution and increasingly negative 
as the orbits become predominantly circular, with the lim- 
iting case of —OQ for purely circular orbits. It is unlikely, 
however, that the full range of allowed values of 13 is covered 
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Figure 6. The projected gravitating mass profile for tlie MMM (dotted line) is compared with that of t he MMMC (solid l ineV Also 



shown is the 9 0% error envelope for the MMMC (dashed lines). Overlaid are the strori g lensing points from Kneib et al. (1995)| and 
&; Mao (1994) (boxes) together with the weak lensing points of ^quires ct al. (1996)[ 



Locb 



after correction for cluster mass (as predicted by the MMMC) within the lensin; 
MMMC provides a good match to both the strong and weak lensing points (especially after re-normalisation of the latter) 



(crosses). Also shown are the weak lensing points 
reference annulus (diamonds). It can be seen that the 




100 
Radius (kpc) 



Figure 7. Derived gas entropy profiles for the inner 500 kpc region are plotted. The MMMC (solid line) can be seen to provide a 
considerably flatter central entropy profile than the MMM model (dash-dot line), where both use a King model for the gas density. If an 
NFW gas density form is used with the MMMC (dashed line), the derived entropy can be seen to increase at all radii. The remaining 
best-fit model profiles (dotted) are included for comparison. 
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in any cluster. On the basis of numerical simulations, Vepe; 
& Dom nguez-Tenreiro (1992) found that over the lifetime of 
letore, anieotropioB corroeponding to a /3 more noga — 



8 SUNYAEV-ZEL'DOVICH EFFECT 



most cl 

five tVi jn -1 R arp nnliVply tn havp Viarl timp tn Hpwlnp 



The Sunyaev-Ze l'dovich microwave decrement (Sunyaev & 



Zol'dovich 1972) results from inverse-Compton scattering of 



For the galaxy surface density profile, a standard 
modified-Hubble profile (E = Eq [l + (fp/rc)^] ), with a 
canonical core radius of 250 kpc, has b een used. The same 
parameterisation has been adopted by Natarajan & Kneib 
(1996)Pwho carried out an analysis rather similar to that 
presented here. The line of sight velocity distribution has 
been fitted by a linear ramp model using a maximum likeli- 
hood 



ffitliod. The best-fit has a central velocit}' dispeisiou 



cosmic background photons by electrons in the cluster gas. 
The magnitude of the efi^ect depends upon the integral of 
the gas pressure along the line of sight and hence has a 
different dependence on gas density than the X-ray surface 
brightness. By combining an analysis of the X-ray emission 
with the observed SZ effect, it is possible to determine the 
distance of the cluster, and hence Hp. Recent mea s urements 
(iBirkinshaw fc Ifughes 1994 [Jones et al. 1993 



of 141?> tf^ km 



yhilo the gradient, although conDiatont 



with zero, is such that the dispersion falls with radius at a 
rate of — 114^240 ^m s~^ arcmin"^. 

The results of the analysis are shown in Fig |^, in which 
the calculated anisotropy parameter profiles for both of 
the most interesting X-ray mass models (the MMM and 
MMMC) are shown. Ahhough the shape of the MMMC 
ani sotropy profile is in detail d ifi'erent from those presented 
by Natarajan & Kneib (1996), the conclusion is the same - 



Saunders 



1996) have been used to place constraints upon the H ubble 
constant in this manner. Birkinshaw & Hughes (1994) mea- 



sure the decrement in a linear strip acr oss t he cluster, deriv- 



ing the ID profile. Jones et al. (1993) and Saunders (1996) 



there is a divergence in the degree of anisotropy in the core 
of A22]| 8. This ia a result of the high central line of sight vc - 



work with 2D images of the decrement, allowing construc- 
tion of a parameterised model for the observed microwave 
decrement. 

Several different observations of the decrement for 
A2218 have been combined with X-ra y data, assuming an 
isothermal plasma, to constrain Hp (Bilk fc White 197^ 



locity d ieporBion being inconcietont with the mace provided 



McHardy et al. 1990t [Birkinshaw fc Hughes 1994 



Jones 



1995). The results obtained vary widely, wit h the most dis- 



by eith er the MMM or MMMC. 



In 



parat e estimates being 241^q km s ^ Mpc ^ ( McHardy et al 



nniaotjopy profiloo, tho aniaotjopy plummoto to 



OQ at bi non aoro radiuc (within tho central " > 1' in tho i 



1994) 



199C) and 65+25 km s"^ Mpc"^ (Birkinshaw fc Hughei 



of the MMMC model) and within this radius the solution 
is unp hycioal) requiring an imaginarjr i i folooity dieporeion i 



In Fig P the predictions of o ur X-ray analysis are 



compa red with the measurements of 



There Jrc turr, pr,agiKlp T-cacr.na fr.r tViia Piratly in rnmpar- 

ison with other clusters, A2218 has an unusual ly large cen- 
tral velocity dispersion. Numerical simulations (Schindler fc 



Birkinshaw fc Hughe; 



Bohrin{;er 1993) have shown that the velocity dispersion can 



increase by up to a factor of 2 during a merger event. This 
occurs during the violent relaxation phase of the merger. At 
this time, the assumptions underlying any analysis based 

tViP Tpana pqnatinn arp rmralrrl Thia prmnrlpa nnp 



(1994) and the allowed envelope of Jones et al. (1993), as- 
suming Ho=50 km s"'^ Mpc"'^. The typical statistical uncer- 
tainty for a single model is similar to the scatter between the 
best-fit model predictions; both of these are small compared 
to the SZ error envelope. 

At radii greater than 100 kpc, the observed and pre- 
dicted decrements are in excellent agreemen t. Howeve r, 



"""ay 



upon 

in whic h the unphysical behaviour of the anisotropv param- 



at smaller radii, the decrement observed by Birkinshaw 
Hughes (1994)1 lies significantly below the X-ray model 



eter can be understood. 

Secondly, the unphysical values of /3 may be indicat- 
ing that the spherically symmetric model and assumptions 
that our solution uses may be in error. This is supported 
by several sets of independent evidence which sugge st that 
the clu ster c ore is disturbed, on the scale of ~ 1'. Kneib 
et al. (1 995) find that the strong le nsing data require a bi- 



predicted profiles. Limitations inherent to beam-switching 
single-dish measurements are that they are prone to base- 
line errors and beam dilution. Thus our an alysis is not based 
upon the 



Birkinshaw fc Hughes (1994) 



modal potential, Markevitch (1997) uses high-resolution X 
ray data to indicate that the cluster strongly deviates from 



data, but instead 
makes use of their results only to allow a comparison with 
earlier studies. 

These problem s are not shared by the observation of 



Jones et al. (1993 



both the 



Birkinshaw & Hughes (1994) 



spherical symmetry in the core and recently Girardi et al 
(1997)|]3howed that the combined galaxy spatial and red- 
shift data indicate the presence of two merging galaxy sub- 
clumps. Thus, there are good reasons to believe that A2218 
has recently undergone a merger event, upsetting the virial 
equilibrium in the cluster core. 

In the case of the MMM model, the solutions are un- 
satisfactory throughout the region for which galaxy data are 
available, such that a more widespread upheaval would be 
required to account for the high velocity dispersion. 



whose calculated e nvelope encompasses 
observations and the 
majority of the predicted profiles for Ho =50 km s"'^ Mpc~^, 
with the MMM being an important exception. This discrep- 
ancy occurs because the MMM requires a steep gas temper- 
ature gradient at small radii (see Fig ^), resulting in a high 
prediction for the central decrement. If the gas is assumed 
to be isothermal (which is not statistically allowed by our 
X-ray data) the predicted SZ decrement at large radius is 
significantly greater than that from the best-fit ASCA mod- 
els (see Fig |9|). However, the isothermal model cannot be 
discriminated against on this basis as it remains consistent 
with the large SZ error envelope. 

If an alternative value of Ho is assumed, the predicted 
decrements obtained from the best-fit X-ray models can be 
varied to achieve consistency with the SZ observations. The 
dependency is such that for a given X-ray surface brightness. 
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Figure 8. The anisotropy parameter for the two mass models discussed in the text is plotted as a function of radius. The anisotropy 
derived using the MMMC mass profile is indicated by the solid line, together with its 90% uncertainty envelope (dashed lines). The 
MMM (dash-dot line) is shown together with its upper error bound (dash-3dot line) only, since the lower bound is unphysical at all radii. 
The MMMC is acceptable at r>l', while the MMM is only allowed at radii beyond the observed edge of the optical data. Both models 
predict predominantly radial orbits at large radius. The /3 = —1.5 limit refers to the constraint from simulations (referred to in the text) 
that anisotropics more negative than —1.5 should not have had time to develop. 
(/) 
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Figure 9. Sunyaev-Zel'dovich decrement profiles predicted from the MMM (solid line) and best-fit cluster r aodels (dotted lines) . together 



with the decr ement which occurs when an isothermal model is used (dashed line). The observed results of Birkinshaw & Hughes (1994) 
(crosses) and [joncs ct al. (1993) (dash-dot la error envelope) are overlaid. A Hubble constant of 50 km s Mpc^-'- is assumed. 
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the predicted SZ decrement varies as Ho° ^ Hence, assumme; 
a higher Hubble constant will lead to a lower X-ray predi cted 



decrement. T he 2D SZ observation of Jones et al. (1993) 



and 



Jones (1995)1 is ideal for such a comparison because, first, it 



avoids the uncertainty inherent in the Birkinshaw & Hughes 
(1994)[]inaly sis and, second, upper and lower bounds to the 
allowed decrement are derived. However, it should be noted 
that these decrements have been analytically parameterised 
to allow a fit to the mosaiced SZ image a nd hence include a 
degree of model dependence ( Jones 1995 ). With this caveat 
in mind, it is possible to determine the range of Ho allowed 
by the SZ and X-ray results. 

Under the requirement that at least one of the ASCA 
best-fit mode ls must be consistent with the results of Jones 
et al. ( [1993)1 , we find that Ho can be limited to the very 
conservative range 37-230 km s~^ Mpc~^. The lower bound 
is obtained by determining what value of Ho is required to 
make the SZ prediction from the LTF model (which predicts 
the lo west central decrem ent of any of the models) equal 
to the Jones et al. (1993' upper bound. The upper bound, 
which is clearly ruled out by other Ho determinations, is ob- 
tained by determining the value of Ho which would decrease 
the SZ prediction from the TTF model (which provides the 
highest centr al decrement of any of the models) such that it 
matches the Jones et al. (1993) lower bound. Using a value 
of Ho outside these bounds results in all of the ASCA best-fit 
models becoming inconsistent with the SZ observations. 

The main conclusion of this analysis is that only weak 
constraints on the Hubble constant can be obtained, with 



even the very low McHardy et al. (1990) result being allowed 
within its errors. 

However, this determination neglects an important ad- 
ditional constraint, namely the gravitating masses extracted 
from lensing analysis. If we consider only the MMM and 
MMMC, which have been shown to be the preferred models 
when strong lensing is taken into account, the allowed range 
of Ho is more tightly constrained. 

The predicted central decrements from the MMM and 
MMMC (for Ho=50km s"^ Mpc"^) are 4.8x10"* and 
respectively. These are compared to the 



Jones et al 



4.0x10 

(1993)Qipper bound of 4.0xl0~*. To achieve consistency be- 
tween these results. Ho must be > 62km s~^ Mpc~^ for the 
MMM and > 50km s"^ Mpc"^ for the MMMC. Lower values 
of the Hubble constai it ensure that the pr edicted decrements 
do not lie within the Jones et al. (1993) bounds. 

If the constraint of isothermality (at the fitted temper- 
ature of 7.9lQ gkeV) is applied, the predicted central decre- 
ment is 3.7xl0~*. When the value of Hq is allowed to vary. 



Mpc' 
km s~^ 



the range allowed by the Jones et al. (1993) 
105 km s"^ i\/r««-i 
Hq = 38lit6 



bounds is 30- 
This is consiste nt with the v alue of 
^ derived by Jones (1995). 



Mpc 



In summary, the constraints currently available from 
SZ observations combined with X-ray measurements are 
too weak to constrain the Hubble constant to greater ac- 
curacy than 37-230 km s~^ Mpc~^. This range accommo- 
dates (within errors) th e previously determin ed values of 
24tfo km s "^ Mpc"^ ([McHardy et al. 1990D , 65^25 km 
s"^ Mpc ~^ (pirkinshaw fc Hughes 1994[ ) and iitf^ km s~^ 
Mpc"^ (Jones 1995). However, when lensing observations 



are introduced (which favour the MMMC model) the Hubble 
constant is required to be greater than 50 km s~^ Mpc~^. 
This result highlights the importance of using both non- 



isothermal gas models and an approach which incorporates 
constraints in addition to those provided by the X-ray and 
SZ observations alone. 



9 DISCUSSION 

We have analysed ROSAT PSPC and ASCA CIS X-ray 
data, fitting spherically symmetric emission models to al- 
low the extraction of cluster properties. The use of ASCA 
data allows the gas temperature structure to be discerned, 
a significant advance over earlier instruments. The analy- 
sis procedure utilises this information by fitting a variety 
of parametric forms for gas density and gas temperature or 
gravitating mass. This avoids restricting the unknown tem- 
perature profile to a single form. 

Since the analysis presented here removes the constraint 
of isothermality, it is important to understand the effect 
that this assumption has. When isothermality is applied to 
Equation ^, it degenerates to: 



Mgrav{r) = 



G/imp 



din pgas{r) 



dlnr 



(10) 



Hence the shape of the mass profile is constrained by the 
gas density profile alone. Since the latter is commonly taken 
to follow a King model, the gravitating mass distribution is 
forced to take the form of an isothermal sphere, with p(r) oc 
r at large radius and flattening within a region determined 
by the gas core radius. 

The main results are: 

1) Comparison with strong lensing indicates that the 
previously reported discrepancy is not present when the 
MMM is used. A consequence of the extra mass required 
within the critical radius is a high central gas temperature, 
which may be related to the disturbed nature of the core. 
So, by relaxing the assumption of isothermality and using 
an appropriate parameterisation for the dark matter distri- 
bution (one which includes a mass cusp), the observed X-ray 
and strong lensing data ca n become consi stent. 

2) The suggestion, by Makino (1996), that the cD con- 
tributes a significant fraction of the cluster mass at small 
radius has been tested. The maximum central mass which 
can be added to the MMM, before the Cash statistic shifts 
beyond the 95% error level, is 1.7x10^ -^ Mp^ (within l OOkpc). 
Thus the massive galaxy required by Makino (1996) is ruled 



out, if the underlying DM distribution follows the form of 
the MMM model. The effect of adding the maximum al- 
lowed central galaxy mass is to extend the DM distribution. 
Because of this, the MMMC is found to be more consistent 
than the MMM with the ou ter two data points extracted 
from the Kncib et al. (1995) strong lensing analysis. 

3) At larger radii, 200-1000 kpc, the projected mass 
profiles of the MMM and MMMC are con sistent with the 
weak lensing results of Squires et al. (1996) , When the slope 



of the mass profile is considered, the MMMC becomes the 
favoured model, as the MMM distribution is considerably 
flatter than the trend of the weak lensing data. It is im- 
portant to recognise that the weak lensing analysis provides 
only a lower limit to the actual mass distribution (since clus- 
ter mass is known to reside in the weak lensing control annu- 
lus). Hence, although the MMMC provides more mass than 
the weak lensing results currently allow, this discrepancy 
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may be resolved when statistical lensing of the background 
galaxy population is observed to greater radii (so that a 
control annulus beyond the edge of the cluster can be used) . 

4) The steep central rise in temperature required by the 
MMM leads to a corresponding increase in entropy. Under 
these conditions the gas is likely to be convectively unsta- 
ble within ~ 120kpc. This is either a reflection of the true 
physical state of the ICM or a model-dependent effect, re- 
lated to the parametric forms assumed for the gas density 
and temperature. Evidence for the former is provided by 
Markevitch (1997), whose HRI analysis indicates the exis- 
tence of central substructure, perhaps as a result of merging 
activity. If this is the case, the gas is likely to have been 
violently shock-heated, such that the assumption of equi- 
librium within the cluster core is no longer secure. On the 
other hand, the problem of convective instability is lessened 
with the MMMC, due to the lower fitted central tempera- 
ture. In addition, replacing the King gas density profile with 
an NFW description can resolve the problem entirely. This 
is a consequence of allowing the density to rise, rather than 
flatten, in the cluster core. 

5) Analysis of galaxy orbits provides another probe of 
the core region of A2218. Using the derived anisotropy of 
the galaxy velocities, it is possible to extract information 
about the gravitational potential and possible presence of 
substructure. Beyond the central 1' region, the velocity dis- 
persion data are consistent with the MMMC. At smaller 
radii, a physically reasonable solution for the anisotropy pa- 
rameter cannot be attained with any of the derived X-ray 
mass distributions. This suggests that the central velocity 
dispersion has been raised to its high observed value dur- 
ing the violent relaxation phase of a merger event. Thus, 
even though the galaxy data are too poorly sampled spa- 
tially to indicate the presence of su bstructure, the velo city 
information suppor ts the analyses of Kneib et al. (1995) and 



Markevitch (1997) 



6) Comparison with 2D Sunyaev-Zel'dovich observa- 
tions indicates that the Hubble constant is only weakly con- 
strained, to a range of 37-230 km s"^ Mpc"^ This is due 
to the high level of uncertainty which occurs when errors 
in the SZ and X-ray data are combined. More restrictive 
constraints on the value of Ho can be obtained if additional 
information is used, such as lensing observations. When the 
MMM and MMMC are combined with the SZ data, low 
values of the Hubble constant are ruled out and we find 
that Ho > 50 km s~^ Mpc~^. This result conflicts with 
t he previously deterrn ined value of 24^^q km s~^ Mpc"'^ 
(McHardy et al. 199C), b ut agrees with the later est imate 



of 65^" km s-i Mpc"^ ( |Birkinshaw fc Hughes 1994] ). The 
differences between these results, and that derived here, are 
dominated by several factors. The first, and most important, 
of these is that an unwarranted assumption has been made 
about the gas temperature in these earlier studies - that it is 
isothermal. This biases the analysis and also results in over 



opti mistic error estimates. Secondl y, McHardy et al. (1990) 
and Birkinshaw & Hughes (1994) both use (different) ID 



SZ measurements, which are prone to baseline uncertain- 
ties. Third, these studies use the less accurate gas density 
information from Einstein, rather than ROSAT. 

Note, however, that these results are all based upon the 
commonly made assumption that the X-ray and SZ obser- 
vations can be made consistent purely by manipulating the 



value of Ho. If the cluster gas, at small radius, is not in 
hydrostatic equilibrium, this methodology may be in error. 



10 CONCLUSIONS 

Combining the above results, it appears that A2218 con- 
sists of two distinct regions. At small ra dii (< 1'), the X - 
ray morphology is extremely disturbed (Markevitch 1997 ) 
with both the strong len sing (Kneib et al. 1995) and galaxy 
data (Girardi et al. 1997) indicating the presence of bimodal 
structure. An X-ray model (the MMMC) consistent with 
the strong lensing data requires the gas temperature to rise 
steeply within the core, possibly at such a high rate that the 
gas is convectively unstable. When this model is combined 
with the observed galaxy data, a physical solution cannot be 
obtained for galaxy orbits within the core. Taken together, 
these results suggest that A2218 has recently undergone a 
merger, shock-heating the gas and disturbing the equilib- 
rium of the components within the cluster. This explains 
the lack of a cooling flow in the system. 

Outside the core, our results show that the data can 
be consistently interpreted on the basis of a cluster in equi- 
librium. The gas temperature falls to < 10 keV, consistent 
with the observed galaxy velocity dispersion, and the gas 
entropy increases with radius. The MMMC provides a mass 
profile consistent with both the outer two strong lensing 
points and the weak lensing mass profile (with the proviso 
that the latter is likely to represent an underestimate of the 
true mass profile). When the MMMC is combined with the 
galaxy velocity data, a physical solution for the galaxy or- 
bits is recovered. The SZ data, which are sensitive to the 
gas pressure at large radius, are consistent with the MMMC 
when Ho >50 km s"^ Mpc"\ 

Thus, while it is premature to regard a model such as 
the MMMC as a complete description of the physical struc- 
ture of A2218, it does appear to explain all the data available 
beyond the disturbed core. To probe further, more detailed 
X-ray and weak lensing observations are required. 

One of the principal conclusions to be drawn from the 
present analysis is that it is dangerous to assume an isother- 
mal ICM without supporting evidence. This can lead to bi- 
ased conclusions and underestimated errors. 
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